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Abstract 

Vanadium dioxide (VO2) undergoes a metal-insulator transition (MIT) at 340 K with the struc- 
tural change between tetragonal and monoclinic crystals as the temperature is lowered. The con- 
ductivity a drops at MIT by four orders of magnitude. The low-temperature monoclinic phase is 
known to have a lower ground-state energy. The existence of a fc-vector k is prerequisite for the 
conduction since the k appears in the semiclassical equation of motion for the conduction electron 
(wave packet). Each wave packet is, by assumption, composed of the plane waves proceeding in 
the k direction perpendicular to the plane. The tetragonal (V02)3 unit cells are periodic along 
the crystal's x-, y-, and z-axes, and hence there are three-dimensional A:-vectors. The period- 
icity using the non-orthogonal bases does not legitimize the electron dynamics in solids. There 
are one-dimensional k along the c-axis for a monoclinic crystal. We believe this decrease in the 
dimensionality of the A;- vectors is the cause of the conductivity drop. Triclinic and trigonal (rhom- 
bohedral) crystals have no fc-vectors, and hence they must be insulators. The majority carriers in 
graphite are "electrons" , which is shown by using an orthogonal unit cell for the hexagonal lattice. 

PACS numbers: 71.30.+h, 81.07.Bc, 73.63.-b 
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Introduction - In 1959 Morin reported his discovery of a metal-insulator transition (MIT) 
in vanadium dioxide (VO2) [1]. Compound VO2 forms a monoclinic (MCL) crystal on the 
low temperature side and a tetragonal (TET) crystal on the high temperature side. When 
heated, VO2 undergoes an insulator-to-metal transition around 340 K, with the resistance 
drop by four orders of magnitude. The phase change carries a hysteresis just as a ferro- 
to-paramagnetic phase change. The origin of the phase transition has been attributed by 
some authors to Peierls instability driven by strong electron-phonon interaction [2] , and by 
other authors to Coulomb repulsion and electron localization due to the electron-electron 



interaction on a Mott-Hubbard picture 



St 



5]. 



A simpler view on the MIT is presented in this Letter. The MLC phase is known to have 
a lower ground-state energy than the TET phase. The existence of /c-vectors is prerequisite 
for the electrical conduction since the k appear in the semiclassical equation of motion for 
the conduction electron (wave packet). The TET (V0 2 )3 unit cells are periodic along the 
crystal's x-, y-, and z-axes, and hence there are three-dimensional (3D) /c-vectors. There 
are ID fc-vectors along the c-axis for a MCL crystal. We show that the MIT occurs since the 
dimensionality of the A;-vectors is reduced from three (3) to one (1) in going from the TET to 
the MCL crystals. We also show that the majority carriers in graphite are "electrons" based 
on the orthogonal unit cells for the hexagonal crystals. In this letter conduction electrons are 
denoted by quotation marked "electrons" ("holes") whereas generic electrons are denoted 
without quotation marks. 

Theory - Following Ashcroft and Mermin 6j, we adopt the semiclassical model of electron 
dynamics in solids. It is necessary to introduce fc-vectors: 

k — k x G x -\- ky&y -\- k z G z , (1) 

where e x , e y , e z are the orthonormal unit vectors, since the A;- vectors are involved in the 
semiclassical equation of motion: 

dk 

hk = h— = q(E + vxB), (2) 

where q is the charge of a conduction electron, and E and B are the electric and magnetic 
fields, respectively. The vector 
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is the electron velocity, where e = e(k) is the energy. 
If we introduce the mass tensor Ai defined by 



l^'l^aSv (4) 



then the equations of motion can be written as 
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The mass tensor M. is symmetric: 

rriij = rriji , (6) 

and can be characterized by the effective masses \m*\. If we choose a Cartesian coordinate 

M 

system along the principal axes of the mass tensor, we can write Eq. ((5]) as [7] 

(\u ■ 

m^ = q{E + vxB) i . (7) 

In this form the Newtonian character of the equations of motion is transparent. If an electron 
is in a continuous energy range (energy band), then it will be accelerated by the electric 
force qE following Eq. ([7]), and the material is a conductor. If the electron's energy is 
discrete and is in a forbidden band (energy gap), it does not move under a small electric 
force, and the material is insulator. If the acceleration occurs only for a mean free time 
(the inverse of a scattering frequency) r, the conductivity a for a simple metal is given by 



Drude's formula 



6]: 



a = q 2 nr/m* , (8) 



where n is the electron density and m* the effective mass. 

For some crystals such as simple cubic (SC), face-centered cubic (FCC), body-centered- 
cubic (BCC), tetragonal (TET) and orthorhombic (ORC) crystals, the choice of the or- 
thogonal (x, y, z)-axes and the unit cells are obvious. The 2D crystals can also be treated 
similarly, only the ^-component being dropped. 

We assume that the wave packet is composed of superposable plane-waves characterized 
by the fc-vectors. The superposability is the basic property of the Schrodingier wave func- 
tion in free space. A MCL crystal can be generated from an ORC crystal by distorting 
the rectangular faces perpendicular to the c-axis into parallelograms. Material plane-waves 
proceeding along the c-axis exist since the (x, y) planes containing materials (atoms) are 



periodic in the ^-direction in equilibrium and can execute small oscillations. It has then 
one-dimensional (ID) fc-vectors along the c-axis. In the x-y plane there is an oblique net 
whose corners are occupied by V's for MCL VO2. The position vector R of every V can be 
represented by integers (m, n), if we choose 

Rmn = max + na 2 , (9) 

where ai and a 2 are non-orthogonal base vectors. In the field theoretical formulation the 
field point r is given by 

r = r' + R mn , (10) 

where r' is the point defined within the standard unit cell. Eq. ( TlDl describes the 2D lattice 
periodicity but does not establish fc-space as explained below. 

To see this clearly, we first consider an electron in a simple square (sq) lattice. The 
Schrodingier wave equation is 

ih-^r) = -^V 2 ^(r) + V{r)^{r) , (11) 

where the potential energy V is periodic: 

V(r + R ( $ n ) = V(r) J (12) 
R^ n = ma x + na y = mae x + nae y , (13) 
(a = lattice constant) . 

If we choose a set of Cartesian coordinates (x, y) along the sq lattice, then the Laplacian 
term in Eq. (jlip is given by 

{ 0^ \ 

V 2 iP(x,y)= (_ + _) il>(x,y). (14) 



dx 2 dy 2 

If we choose a periodic square boundary with the side length Na, N =integer, then there 
are 2D Fourier transforms and (2D) fc-vectors. We now go back to the original rhombic 
system. If we choose the x-axis along either a\ or a 2 , then the potential energy field 
V(r) = V{r + R mn ) is periodic in the x-direction but it is aperiodic in the y-direction. We 
can obtain no 2D periodic boundary condition suitable for Fourier transformation. Then, 
there is no 2D fc-space. If we omit the kinetic energy term, then we can still use Eq. (Q 
and obtain the ground state energy (except the zero point energy). The reduction in the 



dimensionality of the k- vectors from 3D to ID is the cause of the conductivity drop in the 
MIT. The MIT proceeds by domains since the insulator (MCL) phase has the lower degrees 
of symmetry. Strictly speaking, the existence of ID /c-vectors allows the MCL material to 
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measured 



have a small conductivity. This happens for VO2, see below. Wu et al. 
the resistance R of individual nanowires ^N x Yi_ x O x with tungsten (W) ranging up to 1.14 
%. The nanowires are grown with the wire axis matching with the c-axis of the high- 
temperature rutile structure. The transition temperature To decreases from 340 K, passing 
the room temperature, to 296 K as the concentration x changes from to 1.14 %. The 
temperature dependence of R for the low temperature phase is semiconductor-like. That is, 
the resistance R decreases with increasing temperature. The behavior can be fit with the 
Arrhenius law: 

a oc R^ 1 cx exp(— e^/k^T) , (15) 

where e a is the activation energy and k-Q the Boltzmann constant. There are seven (7) 
crystal systems, as seen in AM's book They are: cubic (CUB), TET, ORC, MCL, 
rhombohedral (RHL), hexagonal (HEX), and triclinic systems. There are no k- vectors for 
triclinic and RHL systems. Materials forming these crystals must be true insulators. Arsenic 
(As) and Bismuth (Bi) form RHL crystals, and they are insulators. Diamond (C), silicon 
(Si) and germanium (Ge) form diamond (DIA) crystals. A DIA lattice can be decomposed 
into two FCC sublattices, and can therefore be treated similarly to a CUB crystal. A 
sizable number of elements including graphite form HEX crystals. HEX crystals can be 
treated similar to ORC crystals by choosing orthogonal unit cells. See below for the case of 
graphite. 

Graphene - We consider a graphene which forms a 2D honeycomb lattice. The Wigner- 
Seitz (WS) unit cell,a rhombus (yellow lines) shown in Fig. d(a), contains 2 C's. We showed 



in our earlier work 10j that graphene has "electrons" and "holes" based on the rectangular 
unit cell (black solid lines) s 
details can be found in Refs. 

has the charge — e (+e) and a size of the rectangular unit cell, generated above (below) the 



rown in Fig. [U(b). We briefly review our calculations. More 
llo| . We assume that the "electron" ( "hole" ) wave packet 
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Fermi energy e^. We could show [10J that (a) the "electron" and the "hole" have different 
charge distributions and different effective masses, (b) that the "electrons" and "holes" move 
in different easy channels, (c) that the "electrons" and "holes" are thermally excited with 
different activation energies, and (d) that the "electron" activation energy e\ is smaller than 




(b) 
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FIG. 1. (a) WS unit cell, rhombus (yellow solid lines) for graphene. (b) The orthogonal unit cell, 
rectangle (black solid lines), (c) An orthogonal unit cell (white solid lines) viewed from the top. 
The carbons (circles) in the A (B) planes are shown in blue (orange). 



the "hole" activation energy e 2 ; 



si < e 2 



(16) 



Thus, "electrons" are the majority carriers in graphene. The thermally activated electron 
densities are then given by 



rij(T) = rije £ i/ k & T ; rij = constant, 



(17) 



where j = 1 and 2 represent the "electron" and "hole" , respectively. Magnetotransport ex- 
periments by Zhang et al. ll| indicate that the "electrons" are majority carriers in graphene. 
Thus, our theory is in agreement with experiments. 

Graphite - Graphite is composed of graphene layers stacked in the manner ABAB- • • 
along the c-axis. We may choose an orthogonal unit cell shown in Fig. [T](c). The unit 
cell contains 16 C's. The two rectangles (white solid lines) are stacked vertically with the 



interlayer separation, cq = 3.35 A much greater than the nearest neighbor distance between 
two C's, a = 1.42 A: 

c >a . (18) 

The unit cell has three side-lengths: 

&i = 3a , b 2 = 2^00, h = 2c . (19) 

Clearly, the system is periodic along the orthogonal directions with the three periods 
(b\, b 2 , 63) given in Eq. f lT9|) . We assume that both "electron" and "hole" have the same unit 
cell size. In summary the system is orthorhombic with the sides (pi, b 2 , 03), b\ 7^ b 2 , bi 7^ 63, 
b 2 ^ h- 

The negatively charged "electron" (with the charge — e) in graphite are welcomed by the 
positively charged C + when moving vertically up or downwards. That is, the easy directions 
for the "electrons" are vertical. The easy directions for the "holes" are horizontal. There 
are no hindering hills for "holes" moving horizontally. Hence, the "electron" in graphite has 
the lower activation energy e than the "hole": E\ < e 2 . Then, "electrons" are the majority 
carriers in graphite. The thermoelectric power (Seebeck coefficient) measurements by Kang 



et al. 12] show that the majority carriers in graphite are "electrons", which is in agreement 



with our theory. 

It is often said |l3j| that since the separation distance Co is much greater than the nearest 
neighbor distance do, the conduction in graphite is two-dimensional, and can be discussed 
in terms of the motion in the graphene as a first approximation. We take a different point 
of view. The conduction electrons move as wave packets having the 3D orthogonal unit cell 
sizes. The conduction is two-dimensional because of the inequality (I18jl . But the transport 
behaviors in graphite and graphene are very different because of the different unit cells. 

The construction of the orthogonal unit cell developed here can be followed in other 
material forming HEX crystals. Zinc (Zn) and Beryllium (Be) form HEX crystals. The 
closed orbits on the coronet-like Fermi surface j(J] generate cyclotron resonance, which may 
be discussed using the orthogonal unit cells. 

In summary we established that (a) The MIT in VO2 directly arises from the lattice 
structure change between the TET and the MCL crystals. The TET (MCL) crystal has 3D 
(ID) fc-vectors. The reduction in the dimensionality of the fc-vectors causes a conductivity 
drop, (b) Triclinic and trigonal crystals have no fc-vectors, and hence they are insulators, 



(c) Hex graphite has "electrons" as the majority carriers. 
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